纤维预制件渗透率的预测及其浸润过程有限元模拟

【摘要】先进树脂基复合材料因其具有高比强度、高比模量、耐腐蚀、耐疲劳、可设计性强和可大面积整体成型等诸多优点而成为航空航天领域不可或缺的重要战略材料。然而,复合材料制造高成本限制了其进一步发展,为此低成本的树脂传递模塑(ResinTransferMoulding,RTM)成型已成为航空工业的主流技术之一,但在其制品中经常会出现干斑、孔隙等浸润缺陷,这些缺陷会极大影响产品的质量和性能,而预制件渗透率的不均匀性是引起空隙或干斑的一个重要原因。纤维预制件的渗透率描述了预制件对树脂流动的阻碍程度,是RTM成型过程中树脂流动的一个重要参数。渗透率作为预制件的关键性能之一,与其结构密切相关,由于它与预制件多层次结构的复杂相关性而使得渗透率的预测一直面临着挑战；另外,树脂在纤维中的充分浸润流动是RTM复合材料制品性能的关键影响因素。因此,探讨预制件渗透率与其结构参数的相关性并对其精确预报,进而通过数值模拟预测RTM浸润工艺过程是目前亟待解决的关键问题。本文建立非弯折织物(Non-CrimpFabric,NCF)预制件无缝线和有缝线两类单胞,针对纤维预制件具有双级多孔的特性,在基于双尺度多孔介质的束内／束间耦合流动基础上,数值模拟树脂在预制件细微观结构中的浸润耦合流动行为。采用Brinkman方程描述树脂在纤维束内孔隙的流动,采用Stokes方程描述树脂在纤维束间区域的流动,建立在预制件束内和束问树脂耦合流动的数学模型,利用体积平均化理论,求解速度场和压力场,通过Darcy定律开展NCF预制件面内渗透率的预测。探究了影响渗透率的关键结构参数,开展了不同纤维铺层方式对NCF预制件面内渗透率的影响研究,揭示了纤维预制件的多层次结构参数对面内渗透率的影响规律,在渗透率的预测方面取得了较大进展。对比了有缝线和无缝线两种类型NCF单胞的渗透率大小,重点分析了纤维束形状大小(包括纤维束宽度B、高度h和纤维束沿纤维方向截面椭圆长半轴c、束间距b)及纤维布不同铺层方式等结构因素对预制件面内渗透率的影响。通过引入Morris全局灵敏度分析方法量化了纤维预制件结构参数对渗透率的影响程度,为纤维预制件结构进一步优化提供了重要的理论依据。研究了纤维束在缝制分叉时等效直径、缝线的尺寸和分布、缝线偏离中心距离等因素对纤维预制件面内渗透率的影响。结果表明：没有缝线的预制件单胞面内渗透率大于有缝线的预制件渗透率。预制件渗透率随着纤维束宽度的增加而降低,随着纤维束的高度、纤维束截面边沿椭圆弧的长半轴和纤维束间距的提高而提高。当纤维束宽增加1.2倍时,无缝线的单胞渗透率降低了约52%,而有缝线的单胞渗透率降低了约54%；高度提高一倍时,对应无缝线单胞渗透率增加了约2倍,而有缝线单胞渗透率增加了1.73倍；长半轴增加3.2倍,对应无缝线单胞渗透率增加了30%,有缝线单胞渗透率增加了29%；纤维束间距增加4.7倍,对应无缝线单胞渗透率增加了11.3倍,而有缝线单胞渗透率增加了14倍。预制件渗透率对于不同参数的灵敏度顺序是：b>h>B>c。预制件渗透率对纤维束间距是最为敏感的,对于纤维束截面边沿椭圆弧的长半轴相对不敏感。NCF预制件因缝制时纤维束的分叉阻碍了流体的渗流速度,当流道内交叉束的等效直径与纤维束间距离比是0.7时,纤维束分叉情况下渗透率比没有分叉的预制件单胞渗透率降低约80%。预制件渗透率随着缝线中心位置逐渐偏离流道中心略有增加,随着缝线倾斜角度增加而略有降低。基于预制件纤维束间流道的压力降和其几何结构参数的关系,建立了纤维预制件面内渗透率与其结构相关联的解析模型；对于缝制的单胞,通过引入修正因子和缝线等效方法,建立了纤维预制件渗透率与关键参数的结构相关性解析模型。应用该模型预测的渗透率与有限元计算结果进行对比,结果表明所提出的解析模型在一定范围内能够精确预测纤维预制件面内渗透率,从而实现了多层次结构预制件渗透率的快速精确预测,为完善RTM成型充填理论、优化工艺参数和降低成本起到重要作用。针对北京航空材料研究院先进复合材料重点实验室提出的“离位”增韧技术,揭示了增韧层对预制件渗透率的影响规律和机理,重点分析了增韧层的厚度和渗透率以及不同纤维铺层角度对纤维预制件Z向渗透率的影响。结果表明：纤维预制件Z向渗透率随着增韧层厚度的增加而降低,随着增韧层渗透率的增加而增加；[0]2铺层的纤维预制件Z向等效渗透率远大于其他铺层角度的预制件渗透率,在[0/30]、[0/45]、[0/60]和[0/90]铺层中,[0/45]铺层预制件渗透率最大,[0/30]铺层预制件渗透率最小。纤维体积分数、流道比表面积和流道分布结构决定着渗透率的大小；所提出的数值模拟模型能够精确预测Z向渗透率,作为Z向流动RTM成型模拟的关键参数,它的精确预测对进一步RTM成型的设计和优化提供了重要指导。针对树脂在多尺度纤维预制件中的流动特点,引入了不饱和因子修正Darcy定律,建立了描述树脂在纤维预制件中非稳态流动的偏微分方程(PartialDifferentialEquation,PDE),利用COMSOL软件作为求解器求解RTM成型时树脂流动的控制方程,研究了树脂流动前沿演化规律,将其与解析解和实验结果进行对比,结果表明该模型能够精确预测树脂在纤维预制件中的非稳态流动。针对“离位”增韧技术和Z-RTM成型技术,研究了RTM非稳态浸润过程注射口压力与时间关系的曲线,并与实验结果进行了对比。模拟了树脂在层间未增韧和“离位”增韧纤维预制件束内和束间的流动,实现了树脂在纤维预制件细微观层次浸润的可视化。这种可视化结果为预测树脂在预制件中的宏观流动提供了重要补充并为实际工艺提供了一定指导作用。通过分析发现增韧层使得流动前锋更加平滑,从而减少了成型制品的缺陷。研究了不同成型工艺参数(包括树脂黏度、注射压力和不同渗透率等)对流动前沿的影响规律,从而为实际生产提供了一定的指导作用。
【作者】董抒华；
【导师】王成国；
【作者基本信息】山东大学，材料学，2014，博士
【关键词】树脂传递模塑成型；渗透率；有限元方法；“离位”增韧；非稳态浸润；

【参考文献】
[1]潘丹.中国内地都市报公益活动新闻专版研究[D].南昌大学,新闻传播学（专业学位）,2014,硕士.
[2]张晶.法国电信战略转型:循序渐进式革命[J].世界电信.2006(12)
[3]胡裕效.正常人和2型糖尿病人脑葡萄糖代谢局域及网络属性的初步研究[D].南京大学,内科学,2014,博士.
[4]谢静静.小麦全产业链环节的成本收益研究[D].西南大学,农业经济管理,2014,硕士.
[5]袁万.基于2.5维有限元饱和土地基中隔振沟屏障性能研究[D].浙江大学,岩土工程,2013,硕士.
[6]陈思莹.不锈钢板激光焊接变形的研究[D].大连海事大学,机械工程,2013,硕士.
[7]刘岚.网络环境对高校德育管理的影响及对策研究[D].上海师范大学,教育经济与管理,2004,硕士.
[8]潘伟,井元伟.基于遗传算法的混合降阶H_2/H_∞控制器[J].控制与决策,2005,09:967-970+980.
[9]卢腾超.海洋立管涡激振动的数值模拟研究[D].哈尔滨工业大学,船舶与海洋工程,2014,硕士.
[10]高石磊.热休克蛋白90在急性髓系白血病患者血清中的表达及临床意义[D].河北医科大学,内科学,2014,硕士.
[11]王银河,胡钧,崔乐远,章云.基于参数辨识和T-S模糊模型的一类非线性系统镇定控制器设计[J].控制与决策,2012,06:919-922+928.
[12]谭绍军.HX_N5型内燃机车动力学性能分析[D].西南交通大学,车辆工程,2013,硕士.
[13]李文强.金属氮化物团簇吸附H的密度泛函理论研究[D].兰州理工大学,凝聚态物理,2014,硕士.
[14]越鲜梅.基于图像识别的向日葵叶部病害诊断技术研究[D].内蒙古工业大学,控制理论与控制工程,2013,硕士.
[15]曾建平,程鹏.与“观测器型H_∞控制器设计的LMI方法”一文商榷[J].自动化学报,2001,05:735-736.
[16]王莉.基于改进模糊认知图的入侵检测方法研究[D].天津理工大学,计算机应用技术,2013,硕士.
[17]杨文军.融资融券对股票市场质量影响的实证研究[D].西南交通大学,金融学,2014,硕士.
[18]姚小涛.高校数字图书馆信息服务平台的构建[J].情报杂志,2004,10:106-107.
[19]叶芸.欧债危机对中国出口导向型企业影响的结构性分析及启示[D].华东师范大学,金融学,2013,硕士.
[20]何东轲.有限公司股东控制权滥用问题研究[D].兰州大学,公司法（专业学位）,2013,硕士.
[21]郭万林,张助华,陆鹏.Sp2碳材料中的磁电效应[A].中国材料研究学会.2011中国材料研讨会论文摘要集[C].中国材料研究学会:,2011:2.
[22]罗丹宇.考试焦虑治疗过程中条件性情绪与认知因素的技术干预[D].辽宁师范大学,心理健康教育（专业学位）,2012,硕士.
[23]薛仁经.电子器件高温稳态寿命试验问题研究[J].质量与可靠性.2002(01)
[24]张清.江西省土地整治战略研究[D].江西农业大学,土地资源管理,2012,硕士.
[25]樊宏涛.一类求解非Hermitian鞍点问题的广义局部HSS迭代法[D].兰州大学,计算数学,2014,硕士.
[26]王元刚.面向Agent个体的信念修正系统设计与实现[D].大连海事大学,计算机技术,2014,硕士.
[27]易金容.大鼠高位脊髓损伤后心肌细胞的能量代谢变化和过氧化物体增殖物激活型受体α的表达[D].福建医科大学,麻醉学,2014,硕士.
[28]陈联林.弹性波测试在岩土工程勘察中的应用[J].石油地球物理勘探,1996,S2:27-30+151.
[29]魏晓帅.古钞《文选集注》残存永明诗文研究[D].郑州大学,中国古代文学,2013,硕士.
[30]王跃军.发动机制造企业生产物流优化研究[D].山东大学,工业工程（专业学位）,2012,硕士.
[31]赵恒.甲萘威在甘蓝上的残留分析、膳食风险评估及储藏稳定性的研究[D].河北师范大学,分析化学,2012,硕士.
[32]李慧.基于视频的考场异常行为识别研究[D].天津财经大学,计算机应用技术,2012,硕士.
[33]王江华.东安动力融资收购ABC公司方案设计[D].吉林大学,工业工程,2013,硕士.
[34]王瑗钟.几种锂电池纳米电极材料的制备与电化学性能研究[D].中北大学,化学工程,2013,硕士.
[35]安卓,牛得草,文海燕,杨益,张洪荣,傅华.氮素添加对黄土高原典型草原长芒草氮磷重吸收率及C:N:P化学计量特征的影响[J].植物生态学报,2011,08:801-807.
[36]李刚.北阀公司薪酬体系优化分析[D].中国地质大学（北京）,工商管理,2013,硕士.
[37]姚晨.超小超顺磁氧化铁标记大鼠脂肪源间充质干细胞生物学特性的初步实验研究[D].南方医科大学,神经外科学,2012,硕士.
[38]李幸栗,刘文利,马涛.新疆Y区三维资料速度场研究与应用[J].石油地球物理勘探,1997,S1:75-85+193.
[39]陈杰.中式快餐顾客满意因素实证研究[D].西南大学,企业管理,2013,硕士.
[40]赵阳.基于水足迹理论的昆明市水资源可持续利用研究[D].吉林大学,地质工程,2013,硕士.
[41]王利.融合SIP的TD煤炭专网移动性管理研究与实现[D].西安科技大学,通信与信息系统,2013,硕士.
[42]徐霜.“四月影会”摄影作品艺术价值研究[D].重庆大学,摄影艺术研究,2014,硕士.
[43]徐鹏飞.滑动表面微结构仿真设计及试验研究[D].南京航空航天大学,机械设计及理论,2013,硕士.
[44]刘玉文.颈椎生理弧度的恢复对CSM术后症状改善的影响[D].南京医科大学,临床医学,2013,硕士.
[45]魏文园,张锁喜,姚月清.江苏省减轻自然灾害对策初探[J].地震学刊.1993(02)
[46]田家硕.中西方户外广告设计发展的差异性研究[D].沈阳师范大学,美术学,2014,硕士.
[47]张田文,李仲荣.四值逻辑在图象处理中的应用[J].自动化学报,1989,01:45-49.
[48]邓超,郭茂祖.基于Tri-Training和数据剪辑的半监督聚类算法[J].软件学报,2008,03:663-673.
[49]覃柳师.高校教学管理系统的设计与实现[D].厦门大学,软件工程,2014,硕士.
[50]刘玲绯.Bacillus sp.DSL09-60b菌株的选育及其PLA降解酶的研究[D].东北师范大学,微生物学,2012,硕士.